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(54) INTAKE AIR FLOW ESTIMATE DEVICE FOR INTERNAL COMBUSTION ENGINE 

(57)Abstract: 

PROBLEM TO BE SOLVED: To estimate the intake air flow 
introduced in a combustion chamber through an intake passage 
precisely regardless of an engine operation state. 
SOLUTION: The calculation device 30 of an internal combustion 
engine 10 estimates the temperature after the prescribed term 
from the present point of time of the intake air staying on the down 
streamside part 28 than the throttle valve 24 of the intake passage 
20 based on the fluid model regarding the intake air in an intake 
passage 20 and the temperature information at the present point of 
time of the intake air. The calculation device 30 estimates the 
intake air flow introduced in the combustion chamber 29 from the 
downstream part 28 after the prescribed term from the present 
point of time while reflecting the estimation result. 
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* NOTICES * 

iTPO and NCIPI are not responsible for any 
damages caused by the use of this translation. 

1 .This document has been translated by computer. So the translation may not reflect the original precisely. 

2 **** s hows the word which can not be translated. 
3. In the drawings, any words are not translated. 



CLAIMS 



[Claim(s)] 

[Claim 1] In the inhalation air content prediction equipment of the internal combustion engine which adjusts the 
amount of the inhalation air introduced into a combustion chamber by the throttle valve of an inhalation-of-air 
path To the inhalation air in said inhalation-of-air path The related fluid model The temperature information on 
this time of this inhalation air The amount of the inhalation air introduced into said combustion chamber from 
said downstream part in a predetermined period from this time, presuming the temperature after a 
predetermined period from this time of the inhalation air which is alike, is based and exists in the part of the 
downstream rather than said throttle valve of said inhalation-of-air path, and reflecting this presumed result 
Inhalation air content prediction equipment of the internal combustion engine characterized by having a 
prediction means to predict. 

[Claim 2] In the inhalation air content prediction equipment of an internal combustion engine according to 
claim 1 said prediction means The fluid model built by the heat energy conservation law list about the 
inhalation air of said downstream part based on a law of mass conservation, Based on the pressure and 
temperature in this time of inhalation air of said downstream part, the pressure and temperature after a 
predetermined period are presumed from this time of the inhalation air of said downstream part in the amount in 
this time and temperature of the inhalation air which flows into said downstream part, and a list. Inhalation air 
content prediction equipment of the internal combustion engine characterized by being what predicts said 
inhalation air content based on the this pressure presumed and temperature. 

[Claim 3] In the inhalation air content prediction equipment of an internal combustion engine according to 
claim 2 said prediction means mt (0) and temperature for the amount in this time of the inhalation air which 
flows into said downstream part Tha, Pm (0) and temperature for the pressure in this time of the inhalation air 
of said downstream part Tm (0), When the ratio of specific beat of mc (i) and air is set to kappa and the volume 
of R and said downstream part is set [ the amount of the inhalation air introduced into a combustion chamber 
from said downstream part ] to V for a gas constant, The pressure Pm (i) of the inhalation air of the discrete 
formula (1) obtained from said each conservation law, (2) and said inhalation air content mc (i), and said 
downstream part, and relational expression with temperature Tm (i) (3) 
[Equation 1] 

P mW = ^»<' - » + A ' • * ' Y ' ("M' ~ 0 ' T ha - m c {t - i) . T m (i - i)} • -(1) 

P P R 

(0 = -=T" (f - I) + A/ • — • {m t {i - i) - m c u - i>) - - -(2) 

AtiMWcf^m tfwd = Y,(Atoj) 

f { p m <'>, T m (;)): P m (a T m (o \Z o V * X <D m $C 

Inhalation air content prediction equipment of the internal combustion engine characterized by being what 
presumes the pressure Pm (n) and temperature Tm (n) after the predetermined period tfwd from this time of the 
inhalation air of said downstream part by performing processing only whose count n of predetermined it is alike 
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and repeats the based operation. 

[Claim 4] In the inhalation air content prediction equipment of an internal combustion engine according to 
claim 3 said prediction means Perform said a series of repetitive operation for every predetermined interrupt 
period, and it sets to this interrupt timing. Inhalation air content prediction equipment of the internal combustion 
engine characterized by being what computes either [ at least ] the pressure Pm in this time of the inhalation air 
of said downstream part used for the repetitive operation in next interrupt timing (0), or the temperature Tm (0). 



[Translation done.] 
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DETAILED DESCRIPTION 



[Detailed Description of the Invention] 
[0001] 

[Field of the Invention] This invention relates to the inhalation air content prediction equipment of the internal 
combustion engine which predicts the amount of the inhalation air introduced into an engine combustion 
chamber. 
[0002] 

[Description of the Prior Art] In case the amount (inhalation air content) of the air introduced into an internal 
combustion engine's combustion chamber through an inhalation-of-air path determines the controlled variable 
of internal combustion engines, such as fuel-injection control and ignition- timing control, as what reflects 
engine operational status most with an engine rotational frequency, it is used, and this is usually called for based 
on the detecting signal from an air flow meter or an intake-pressure sensor prepared in the inhalation-of-air 
path. 

[0003] By the way, it faces computing an internal combustion engine's controlled variable based on an 
inhalation air content, and if the inhalation air content in this time called for from detecting signals, such as an 
air flow meter, is used as it is, and it is at the time of transient operation, this controlled variable will come to be 
determined based on a different inhalation air content from the event of the controlled variable being actually 
reflected in engine control. Thus, if the inhalation air content used for the decision of a controlled variable 
differs from the inhalation air content in case the controlled variable is reflected in engine control, control 
adapted to actual engine operational status cannot be performed, but even if it raises the control precision, a 
limitation will come to be generated naturally. 

[0004] Then, the inhalation air content after a predetermined period is predicted from this time, and the 
equipment which determined the controlled variable based on the forecast is conventionally proposed so that it 
may be indicated by JP,2-42160,A. By performing such prediction, the inhalation air content at the time of a 
controlled variable being reflected in engine control can be foreseen beforehand, this controlled variable can be 
determined now, and improvement in control precision comes to be achieved. 

[0005] Moreover, an inhalation air content changes also with an engine rotational frequency, or intake-air 
temperatures besides the opening of a throttle valve. For example, since the consistency of inhalation air 
becomes large so that an intake-air temperature becomes low, an inhalation air content also comes to increase. 
Therefore, in case an inhalation air content is predicted as mentioned above, it becomes important when 
obtaining a prediction result accurate [ also taking into consideration the effect of such an intake-air 
temperature ]. 

[0006] He detects the intake-air temperature in the location of the upstream, and is trying to make that intake-air 
temperature (for this to be hereafter, called "upstream intake-air temperature", when it is necessary to 
distinguish especially from the intake-air temperature in other locations of an inhalation-of-air path) reflect in 
prediction of an inhalation air content rather than the throttle valve of an inhalation-of-air path with this point 
and the above-mentioned equipment. 
[0007] 

[Problem(s) to be Solved by the Invention] However, although the above-mentioned conventional equipment is 
making the intake-air temperature reflect in prediction of an inhalation air content, it cannot but perform the 
prediction to the last by making for location-change and its time change of an intake-air temperature not to arise 
during a prediction period in an inhalation-of-air path into a prerequisite. 
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[0008] If surely it was at the time of steady operation by which the opening of a throttle valve is held uniformly, 
since it was in the range which most location-change of the intake-air temperature in such an inhalation-of-air 
path does not have, and can also disregard change of a time intake-air temperature, though the inhalation air 
content was predicted under the above-mentioned prerequisite (i.e., ****** [ only based on the upstream intake- 
air temperature mentioned above, it predicts an inhalation air content ]), gross errors do not arise in the 
prediction result. 

[0009] However, if it is when the opening of a throttle valve changes rapidly at the time of transient operation, 
intake-air temperatures come to differ according to the location in an inhalation-of-air path, and the time change 
cannot be disregarded greatly, either. Therefore, the above-mentioned prerequisite does not meet the actual 
temperature condition of the inhalation air in such an inhalation-of-air path. 

[0010] For example, when the amount of the inhalation air which flows into the downstream part of this valve 
when the opening of a throttle valve increases increases rapidly, the inhalation air of this downstream part will 
be compressed and the temperature (henceforth a "downstream intake-air temperature") will come to rise more 
nearly temporarily than the upstream intake-air temperature mentioned above. Since it will be conversely 
expanded by the inhalation air of this downstream part when the amount of the inhalation air which flows into 
said downstream part on the other hand when the opening of a throttle valve decreases decreases rapidly, a 
downstream intake-air temperature comes to fall more nearly temporarily than an upstream intake-air 
temperature. 

[001 1] thus, temporary, if it is at the time of transient operation from which the opening of a throttle valve 
changes — an imitation — a downstream intake-air temperature can become temperature which is changed 
sharply in time and is different from an upstream intake-air temperature. Consequently, if it is in the 
conventional equipment which predicts by making for change [ be / location-/ it / time ] of an intake-air 
temperature not to arise in an inhalation-of-air path into a prerequisite, so to speak, the effect by the dynamic 
behavior of the intake-air temperature at the time of such transient operation is not reflected in prediction of an 
inhalation air content, therefore the error of the prediction result cannot be disregarded, either. 
[0012] This invention is made in view of such the actual condition, and that object is in offering the inhalation 
air content prediction equipment of the internal combustion engine which can predict the amount of the 
inhalation air introduced into a combustion chamber through an inhalation-of-air path with a sufficient precision 
even if it is at the transient operation time. 
[0013] 

[Means for Solving the Problem] The means and its operation effectiveness for attaining the above-mentioned 
object are indicated below. In the inhalation air content prediction equipment of the internal combustion engine 
which adjusts the amount of the inhalation air with which invention indicated to claim 1 is introduced into a 
combustion chamber by the throttle valve of an inhalation-of-air path To the inhalation air in said inhalation-of- 
air path The related fluid model The temperature information on this time of this inhalation air The amount of 
the inhalation air introduced into said combustion chamber from said downstream part in a predetermined 
period from this time, presuming the temperature after a predetermined period from this time of the inhalation 
air which is alike, is based and exists in the part of the downstream rather than said throttle valve of said 
inhalation-of-air path, and reflecting this presumed result He is trying to have a prediction means to predict. 
[0014] A prediction means presumes the temperature of the inhalation air of said downstream part based on the 
fluid model about the inhalation air in an inhalation-of-air path, and the temperature information (physical 
quantity which has these or correlation, such as a mode of temperature or a temperature change) on this time of 
this inhalation air, and makes the presumed result reflect in prediction of an inhalation air content in the above- 
mentioned configuration. Therefore, even if it is the case where the temperature of this downstream part 
changes from this time transitionally, when that transitional change has been grasped, an inhalation air content 
can be predicted, and a very accurate prediction result can be obtained. 

[0015] Moreover, since the amount of the inhalation air introduced into a combustion chamber through an 
inhalation-of-air path tends to increase, so that the pressure of the inhalation air of the above-mentioned 
downstream part becomes high, and so that the temperature becomes low, it can ask for this based on the above- 
mentioned pressure and temperature. Therefore, the amount of the inhalation air introduced into a combustion 
chamber in a predetermined period from this time based on the pressure and temperature which are they- 
presumed can be predicted now by presuming the pressure and temperature after a predetermined period from 
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this time of the inhalation air of a downstream part. 

[0016] Furthermore, the pressure and temperature after a predetermined period can presume these from this 
time of the inhalation air of this downstream part by relating with physical quantity called the pressure and 
temperature of inhalation air in this time of that downstream part at the amount in this time and temperature of 
the inhalation air which flows into a downstream part at the heat energy conservation law list about the 
inhalation air of a downstream part based on a law of mass conservation, and a list. 

[0017] Invention indicated to claim 1 is set to the above-mentioned inhalation air content prediction equipment 
as a configuration materialized more like invention which indicated to claim 2. Therefore, the - aforementioned 
prediction means The fluid model built by the heat energy conservation law list about the inhalation air of said 
downstream part based on a law of mass conservation, Based on the pressure and temperature in this time of 
inhalation air of said downstream part, the pressure and temperature after a predetermined period are presumed 
from this time of the inhalation air of said downstream part in the amount in this time and temperature of the 
inhalation air which flows into said downstream part, and a list. The configuration that it is what predicts said 
inhalation air content based on the this pressure presumed and temperature is employable. 
[0018] Invention indicated to claim 3 is set to the inhalation air content prediction equipment of an internal 
combustion engine according to claim 2. Moreover, said prediction means mt (0) and temperature for the 
amount in this time of the inhalation air which flows into said downstream part Tha, Pm (0) and temperature for 
the pressure in this time of the inhalation air of said downstream part Tm (0), When the ratio of specific beat of 
mc (i) and air is set to kappa and the volume of R and said downstream part is set [ the amount of the inhalation 
air introduced into a combustion chamber from said downstream part ] to V for a gas constant, The pressure Pm 
(i) of the inhalation air of the discrete formula (1) obtained from said each conservation law, (2) and said 
inhalation air content mc (i), and said downstream part, and relational expression with temperature Tm (i) (3) 
[0019] 
[Equation 2] 

P m (O = P m H - 0 + Ar . K ■ — . {m t 0 - l) • T ha - m c {i - i) - T m (i - i))- - •(]) 
P P R 

-~-(t) = yr-V - i> + Ar - — • (/7f,(> - l) - m c (r - 1)) •••(2) 
*cW = f( P w 0),T m io) "(3) 

f(P m (/), T m in): P m to, T m (/)Ko ^X<0 m & 

By performing processing only whose count n of predetermined it is alike and repeats the based operation, it is 
supposed that it is what presumes the pressure Pm (n) and temperature Tm (n) after the predetermined period 
tfwd from this time of the inhalation air of said downstream part. 

[0020] Since according to the above-mentioned configuration it is discretized and expressed as said fluid model 
is shown in each above-mentioned formula (1) and (2), the pressure Pm (n) and temperature Tm (n) after a 
predetermined period can be presumed based on a comparatively easy algorithm from this time of the inhalation 
air of said downstream part. 

[0021] It sets to the inhalation air content prediction equipment of an internal combustion engine according to 
claim 3 like [ in performing such data processing ] invention indicated to claim 4. Moreover, the - 
aforementioned prediction means Perform said a series of repetitive operation for every predetermined interrupt 
period, and it sets to this interrupt timing. If the configuration that it is what computes either [ at least ] the 
pressure Pm in this time of the inhalation air of said downstream part used for the repetitive operation in next 
interrupt timing (0) or the temperature Tm (0) is adopted Not using a sensor etc., ** can also search for serially 
the pressure Pm of the above-mentioned inhalation air (0), or its temperature Tm (0) through the above- 
mentioned repetitive operation, and simplification of the configuration of inhalation air content prediction 
equipment can be attained now. 

[0022] Moreover, generally, since the responsibility of a general-purpose temperature sensor is low, when the 
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temperature Tm of inhalation air (0) is detected using this sensor, it is difficult to detect change of this 
temperature Tm (0) with a sufficient precision, and it cannot finish being avoided that an error arises in the 
detection result. 

[0023] In this point and the above-mentioned configuration, if especially the temperature Tm of an inhalation 
air content (0) is computed, while being able to attain simplification of a configuration, lowering of the 
predictability resulting from the detection error of such a sensor can be avoided. 
[0024] 

[Embodiment of the Invention] The 1st operation gestalt of this invention is explained below [the 1st operation 
gestalt]. 

[0025] He is trying to apply the inhalation air content prediction equipment concerning this invention to the 
internal combustion engine which has an air flow meter as a sensor which detects an inhalation air content with 
this operation gestalt. First, with reference to drawing 1 , this internal combustion engine's 1 0 outline and the 
configuration of this prediction equipment are explained. 

[0026] An internal combustion engine's 1 0 inhalation-of-air path 20 is constituted by the inlet port 23 grade 
which connects the surge tank 22 connected to the aisleway of a throttle body 21 in which the throttle valve 24 
was formed, and the downstream of this throttle body 21 and this surge tank 22, and a combustion chamber 29. 
After metering of the inhalation air which flows the inhalation-of-air path 20 is carried out by the throttle valve 
24, it is introduced in a combustion chamber 29 at the time of valve opening of an inlet valve 26. 
[0027] Inhalation air content prediction equipment predicts the amount of the inhalation air introduced in a 
combustion chamber 29 in a predetermined period from this time, and is greatly constituted by the various 
sensors 41-45 including an arithmetic unit 30 and an air flow meter 41 . 

[0028] The arithmetic unit 30 is equipped with the operation part 31 which performs data processing, and the 
memory 32 the various function data used on the occasion of the program concerning such data processing or 
its activation were remembered to be. This arithmetic unit 30 is constituted by the electronic control which 
performs various control of an internal combustion engine 10. 

[0029] Others, an intake temperature sensor 42, the throttle sensor 43, the cam angle sensor 44, and crank angle 
sensor 45 grade are contained in sensors 41-45. [ air flow meter / 41 / above-mentioned ] All, the air flow meter 
41 and the intake temperature sensor 42 are formed in the part of the upstream from the throttle valve 24 at the 
inhalation-of-air path 20, and detect the amount and temperature of the inhalation air which flows in a surge 
tank 22 through this valve 24 (henceforth "the upstream intake-air temperature Tha"), respectively. 
[0030] Moreover, with this operation gestalt, the thing of the heat ray type equipped with the detecting element 
410 which consists of a heat ray as the above-mentioned air flow meter 41 is used. Drawing 2 shows the cross- 
section structure of this detecting element 410. As shown in this drawing, this detecting element 410 is 
constituted by the heat ray 412 which consists of platinum etc., and the glass layer 414 which covers the 
perimeter of this heat ray 412. 

[003 1 ] The throttle sensor 43 is formed near the throttle valve 24, and detects the opening (throttle opening TA) 
of this valve 24. The crank angle sensor 45 is formed [ near / where the cam angle sensor 44 carries out closing 
motion actuation of the inlet valve 26 / the cam shaft (graphic display abbreviation) ], respectively near the 
crankshaft (graphic display abbreviation) which carries out revolution actuation of this cam shaft. An arithmetic 
unit 30 incorporates the output signal of these cam angle sensor 44 and the crank angle sensor 45, and detects 
the closing motion stage (valve timing VT) of an inlet valve 26, respectively in the rotational speed (engine 
engine speed NE) and its rotation phase (crank angle CA) list of a crankshaft. In addition, the device (graphic 
display abbreviation) in which the closing motion stage of the above-mentioned inlet valve 26 is changed based 
on engine operational status, such as the engine rotational frequency NE and an engine load, is prepared for the 
internal combustion engine 10 of this operation gestalt. Therefore, the above-mentioned valve timing VT will 
be suitably changed according to engine operational status. 

[0032] Next, prediction processing of the inhalation air content by the equipment of such this operation gestalt 
is explained. The principle which will be the requisite for this prediction processing first is explained. 
[0033] The space from the inhalation-of-air path 20 to a combustion chamber 29 from a throttle valve 24 First, 
the part of the upstream From 27 and this throttle valve 24, rather than an inlet valve 26 by the downstream The 
part of the upstream (It is hereafter called an "upstream part") (It is hereafter called a "downstream part") It 
divides into 28 and three parts called a combustion chamber 29, and the fluid model about the inhalation air 
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which exists in the downstream part 28 among each [ these ] parts 27-29 is built. 

[0034] Namely, the amount of the inhalation air which flows out of "M" and this downstream part 28 the mass 
of the inhalation air which exists in this downstream part 28 in a combustion chamber 29 per time amount If the 
amount (henceforth "throttle-valve through put") of the inhalation air which flows (it is hereafter called "the 
inhalation air content in a cylinder") into the downstream part 28 from "mc" and the upstream part 27 is set to 
"mt" Variation [ per time amount of the above-mentioned mass M ] d (M)/dt can be expressed like a degree type 
(5) based on a law of mass conservation. 
[0035] 
[Equation 3] 

— M = m t - m c " *( 5 ) 
at 

On the other hand, the heat energy in which this inhalation air has them when "Tm", its isochore specific heat, 
and the isotonic specific heat are set to "Cv" and "Cp", respectively serves as "Cv-M-Tm" in the temperature 
(henceforth a "downstream intake-air temperature") of the inhalation air which exists in the downstream part 28. 
Moreover, since the heat energy which the inhalation air which flows into a combustion chamber 29 from "Cp- 
mt-Tha" and the downstream part 28 has [ the heat energy which the inhalation air which flows into the 
downstream part 28 from the upstream part 27 has ] is "Cp-mc-Tm", the variation d (Cv-M-Tm)/dt per time 
amount of above-mentioned heat energy Cv-M-Tm can be expressed like a degree type (6) based on a heat 
energy conservation law. 
[0036] 
[Equation 4] 

^-(CV ' M ♦ T m ) - C P - m t - Tha ~ C p - mc * T m "'(6) 

Moreover, the equation of state about the inhalation air which exists in the downstream part 28 will become like 
the following equation (7), if the pressure (henceforth a "downstream intake pressure") of the inhalation air for 
"V" and the said division 28 is set to "Pm" and it sets a gas constant to "R" for the volume of the downstream 
part 28. 
[0037] 
[Equation 5] 

p m V = Mlt-T m '"(7) 

Here, if the ratio of specific beat of air is set with "kappa" (=Cp/Cv), a formula (5) can express the following 
formulas (9) and a formula (6) like the following formulas (8) using the above-mentioned formula (7). 
[0038] 
[Equation 6] 



Furthermore, if discrete spacing is set to "deltat", the left part of each [ these ] type (8) and (9) can be expressed 
as shown in a degree type (10) and (11). 
[0039] 
[Equation 7] 



Therefore, the following discrete formulas (12) and (13) are drawn from each [ these ] type (8) - (1 1). In 
addition, in this formula (12) and (13), "i" in a parenthesis shows a value and "i-1" shows a value last time this 
time, respectively. 




) •■•(*) 
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[0040] 
[Equation 8] 

F m(0 = O f - i) + - k • — • (m,</ - i) • r Aa - w c (/ - i) ■ T m v - i))- • -(12) 
P P R 

-jr-(0 = -TjP-O - D+ (/w,<f- i) - w c (i - i)) ""(13) 

On the other hand, the inhalation air content mc in a cylinder can be calculated using the following empirical 
formula (14). 
[0041] 
[Equation 9] 

m c = ^~ {a P m -b) .-(14) 

In a top type (14), each of "a" and "b" is constants, and is set as the flow of the inhalation air in the part from the 
downstream part 28 to a combustion chamber 29 at the affecting engine control parameter, i.e., here, based on 
the engine engine speed NE and valve timing VT. The relation between these constants a or a constant b and the 
engine engine speed NE, and valve timing VT is beforehand defined by experiment, and is memorized by the 
memory 32 of an arithmetic unit 30 as a function map. Moreover, "To" is the reference temperature (constant) 
of the inhalation air when setting up each constants a and b by such experiment. 

[0042] As shown in each following formula (15) and (16), if this empirical formula (14) is used, each above- 
mentioned formula (12) and (13) can deform. 
[0043] 

[Equation 10] 

p mM - P m H - i) + A/ • if y- (m t a - i> • 7^ - r 0 • (a ■ P m o - 0 - b)) ---(15) 
Zk (0 = *k (i _ 0 + At.y(m,v - i) - _A_.(a • P m a - i) -6)1 --(16) 

mm > m " " * 

Moreover, it can ask for the downstream intake-air temperature Tm using the formula (17) of the following 

obtained from each [ these ] type (15) and (16). 

[0044] 

[Equation 11] 

T m <n--^- —07) 



Each [ these ] type (1 5) So that clearly from - (1 7) Throttle- valve through put mt in this time, It is based on the 
downstream intake pressure Pm, the downstream intake-air temperature Tm, and the upstream intake-air 
temperature Tha. The downstream intake pressure Pm and the downstream intake-air temperature Tm after the 
above-mentioned discrete spacing deltat can be presumed after a predetermined period (i.e., here) from this 
time. Furthermore, by substituting for a formula (14) the downstream intake pressure Pm presumed in this way 
and the downstream intake-air temperature Tm, the inhalation air content mc in a cylinder after discrete spacing 
deltat can be presumed from this time. 

[0045] Moreover, he is trying to calculate the above-mentioned throttle-valve through put mt in prediction 
processing of this operation gestalt according to two approaches the followings differ. Hereafter, the approach 
at the time of calculating this throttle- valve through put mt is explained. 

[0046] First, the procedure at the time of calculating the throttle- valve through put mt based on the detecting 
signal of an air flow meter 41 is explained. As everyone knows, an air flow meter 41 detects tales doses using 
the heating value taken from this heat ray 412 changing according to an inhalation air content, in case inhalation 
air passes through the perimeter of a heat ray 412. The relation between the detecting signal of this air flow 
meter 41 and an inhalation air content is beforehand called for by experiment, and is memorized by the memory 
32 of an arithmetic unit 30 as a function map. Therefore, an inhalation air content can be calculated based on the 
detecting signal and this function map of an air flow meter 41 . 
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[0047] However, as compared with a heat ray 412, generally, since heat capacity is large, the above-mentioned 
glass layer 414 has the inclination for the temperature Tgs to be overdue to change of the amount of the 
inhalation air which passes near the detecting element 410, and to change. For this reason, if there is an 
inhalation air content calculated from the above-mentioned function map if it is at the time of steady operation 
from which an inhalation air content does not change at the time of transient operation from which a true 
inhalation air content and this inhalation air content of a match change, the inhalation air content calculated 
from a function map comes to shift from a true inhalation air content. 

[0048] Then, in order to foresee beforehand the response delay about the temperature change of such a glass 
layer 414 and to calculate a true inhalation air content, he builds the thermal equilibrium model about this layer 
414, and is trying to compute an inhalation air content with this operation gestalt based on the reverse model. In 
addition, in the following explanation, a true inhalation air content is made into the throttle-valve through put 
mt, and the inhalation air content acquired based on the detecting signal and the above-mentioned function map 
of an air flow meter 41 is distinguished as a map operation value GA. 

[0049] First, if temperature of a heat ray 412 is set to "Tht" (refer to drawing 2 ), the variation d (Tgs)/dt of the 
temperature Tgs of the glass layer 414 can be expressed like the following formulas (18) from the balance with 
the heating value transmitted to the glass layer 414 from a heat ray 412, and the heating value transmitted to 
inhalation air from this layer 414. 
[0050] 

[Equation 12] 

A Tt T * s = M 7 *' - T **)-( c + D J™) T >™) '"< 18 > 

Here, "A", "B", M C", and "D" are constants determined according to **, such as the cross section of a heat ray 
412, die length and its resistivity, and a heat transfer rate between the glass layer 414 and a heat ray 412, a heat 
transfer rate between the glass layer 414 and inhalation air. 

[0051] Moreover, at the time of steady operation, since transfer of the heat between the glass layer 414, and a 
heat ray 412 and inhalation air is lost, as shown in a degree type (19) and (20), the left part d (Tgs)/dt of 
temperature changes of a formula (18), i.e., the amounts of the glass layer 414, is set to "0", and the map 
operation value GA and the throttle- valve through put mt become equal. 
[0052] 

[Equation 13] 
GA = mt —(20) 

Therefore, the following formula (21) is obtained from each above-mentioned formula (19) and (20). 
[0053] 

[Equation 14] 



Tgs - T ha B 



"(21) 



Furthermore, if a formula (1 8) is transformed using this formula (21), the degree type (22) from which the 

temperature Tgs of the glass layer 414 was eliminated will be obtained. 

[0054] 

[Equation 15] 

^-(tM -*=)-(JE-<fiH) 

Moreover, the left part of this formula (22) can be expressed like a degree type (23), if discrete spacing is set to 

"deltat." 

[0055] 

[Equatio n 16] 

d4GA _ JGAW ~ yjGAjj - l) 

http://www4.ipdl.ncipi.go.jp/cgi-bin/tran__web_cgi_ejje 1 0/1 7/2005 



JP ,2001 -041 095,A [DETAILED DESCRIPTION] 



Page 8 of 15 



Therefore, the following discrete formulas (24) are drawn from each above-mentioned formula (22) and (23). In 
addition, in this formula (24), "alpha" and "beta" are constants which become settled by above-mentioned 
constant A-D, and "j" in a parenthesis shows a value and "j-1" shows the value last time this time, respectively. 
[0056] 

[Equation 17] 



The throttle-valve through put mt in this time (j) can be presumed based on the map operation value GA in this 
time (j), and the map operation value GA (j-1) in front of discrete spacing del tat so that clearly from this 
formula (24). 

[0057] Next, the procedure at the time of calculating the throttle- valve through put mt based on engine 
operational status, such as the throttle opening TA, is explained. When a part for opening of a throttle valve 24 
is assumed to be a kind of orifice, the amount (throttle-valve through put mt) of the fluid which passes this 
orifice is expressed by each following formula (25) and (26) as everyone knows. 
[0058] 

[Equation 18] 



Here, "mu" is the flow coefficient of an orifice (a part for opening of a throttle valve 24), and "Pa" is 
atmospheric pressure (constant value). Moreover, "At" is the opening area of a throttle valve 24, and can ask for 
this uniquely from the throttle opening TA. 

[0059] Here, if the time of steady operation by which the throttle opening TA is held uniformly is assumed, a 

formula (25) will become like a degree type (27). 

[0060] 

[Equation 19] 

"a*-M'A t . , \ *( P »TA>*>a) --(27) 

Here, "mtTA" and "PmTA" show the throttle-valve through put mt at the time of steady operation, and the 
downstream intake pressure Pm, respectively. The throttle-valve through put mtTA at the time of these steady 
operations and the downstream intake pressure PmTA can be determined as the flow of the inhalation air in the 
part from the upstream part 27 to a combustion chamber 29 at the affecting engine control parameter, i.e., here, 
based on the throttle opening TA, the engine engine speed NE, and valve timing VT. The relation between the 
throttle-valve through put mtTA at the time of these steady operations or the downstream intake pressure 
PmTA, the throttle opening TA and the engine engine speed NE, and valve timing VT is beforehand called for 
by experiment, and is memorized by the memory 32 of an arithmetic unit 30 as a function map. 
[0061] If this formula (27) is used, a formula (25) can be expressed like a degree type (28). 
[0062] 

[Equation 20] 
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Furthermore, if it is at the time of steady operation, since the throttle-valve through put mtTA and the inhalation 
air content mc in a cylinder become equal (mtTA=mc) and the downstream intake-air temperature Tm and the 
upstream intake-air temperature Tha become equal (Tm=Tha), the above-mentioned formula (14) can be 
expressed like a degree type (29). 
[0063] 

[Equation 21] 

'ha 

And if this formula (29) is used, the above-mentioned formula (28) can be expressed still like a degree type 

(30). 

[0064] 

[Equation 22] 



The throttle- valve through put mt can be presumed based on the value PmTA of the downstream intake pressure 
Pm when assuming that it is at the steady operation time (TA, NE, VT) in the upstream intake-air temperature 
Tha and the downstream intake pressure Pm, the atmospheric pressure Pa, and list in this time so that clearly 
from this formula (30). 

[0065] He is trying to predict the inhalation air content mc in a cylinder after a predetermined period to be the 
throttle-valve through put mt calculated based on each option in this way from this time with this operation 
gestalt based on each formula (14) - (17) mentioned above. In addition, in the following explanation, the 
throttle-valve through put mt asked for the throttle-valve through put mt calculated based on engine operational 
status, such as the throttle opening TA, based on the detecting signal of the 1 air flow meter throttle-valve 
through put estimate mt 41 is distinguished as throttle- valve through put measured value mt2. 
[0066] Next, drawing 3 and drawing 4 are combined, referred to and explained about the detail of the prediction 
processing by the equipment of this operation gestalt. Drawing 3 R> 3 and drawing 4 are flow charts which 
show an example of the procedure at the time of predicting the inhalation air content mc in a cylinder, this 
processing of a series of is performed by the arithmetic unit 30 the interrupt period (for example, — "— every 
8msec. M ) for every predetermined time. 

[0067] On the occasion of this processing, the prediction period tfwd at the time of predicting an inhalation air 
content is set up first (step 100). This prediction period tfwd will be the fuel oil consumption in the last fuel 
injection, and a thing set up according to fuel injection duration, if it puts in another way. Therefore, the 
inhalation air content mc in a cylinder after the time amount which the last fuel injection took from this time 
passes will be predicted here. Incidentally, this fuel injection duration is computed through processing 
concerning fuel injection with this another processing. Moreover, in processing concerning this fuel injection, 
fuel oil consumption is computed based on the newest inhalation air content mc in a cylinder computed through 
this processing. 

[0068] Next, the downstream intake pressure Pm after the prediction period tfwd, the downstream intake-air 
temperatures Tm, and those ratios Pm/Tm are computed based on said throttle- valve through put estimate mtl 
from this time (step 200). In addition, the downstream intake pressure Pm after the prediction period tfwd 
computed below based on the throttle-valve through put estimate mtl the downstream intake-air temperatures 
Tm and those ratios - Pm/Tm with "Pmltfwd", "Tmltfwd", and "Pml tfwd/Tmltfwd", respectively moreover, 
the downstream intake pressure Pm in the middle of the calculation, the downstream intake-air temperatures 
Tm, and those ratios - it writes "Pml (i)", "Tml (i) M , and "Pml/Tml (i)", respectively. [ Pm/Tm ] 
[0069] On the occasion of calculation of this downstream intake-pressure Pmltfwd and downstream intake-air 
temperature Tmltfwd, processing only whose count of predetermined repeats the operation based on each 
above-mentioned formula (15) - (17) and a formula (30) is performed. 

[0070] That is, as shown in drawing 4 , the increment of the counter value i showing the count of a loop of this 
processing (initial value : "0") is carried out first (step 202), and it is judged whether next the prediction periods 
tfwd are del tat 1 or more predetermined time (step 204). Here, the above-mentioned predetermined time deltatl 
is set up equally (for example, "8msec") to the interrupt period of this routine. And when the prediction period 
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tfwd is judged to be these deltatl or more predetermined time (step 204: YES), discrete spacing deltat of each 
above-mentioned formula (15) and (16) is set up equally to predetermined time deltatl (step 206). 
[0071] next, based on a formula (30), the throttle- valve through put estimate mtl (i) computes — having — 
further — after that and each formula (15) - (17) — being based — the downstream intake pressure Pm 1 (i), the 
downstream intake-air temperatures Tml (i), and those ratios — Pml (i)/Tml (i) is computed, respectively (step 
208). And predetermined time deltatl is subtracted from the current prediction period tfwd, and the subtraction 
value (tfwd-delta tl) is set up as a new prediction period tfwd (step 210). 

[0072] Next, when it is judged whether this new prediction period tfwd is below "0" (step 212) and this 
prediction period tfwd is judged to be longer than "0" (step 212: NO), processing after step 202 is performed 
again. Under the present circumstances, if it is judged that the prediction period tfwd is less than deltatl 
predetermined time in step 204 (step 204: NO), discrete spacing deltat will be set up equally to the prediction 
period tfwd (step 207), and processing after step 208 will be performed. 

[0073] And by repeating processing of steps 202-212 If the prediction period tfwd becomes below "0" (step 
212: YES) the downstream intake pressure Pm 1 (i), the downstream intake-air temperature Tml (i), and a ratio 
— Pml /Tml (i) of them downstream intake-pressure Pml tfwd after the prediction period tfwd, downstream 
intake-air temperature Tmltfwd, and a ratio ~ it is set up as Pml tfwd/Tmltfwd, respectively (step 214). 
[0074] For example, when the interrupt period of this routine is " 8msec." and the prediction period tfwd is set as 
"30msec. 11 through processing of previous step 100, Discrete spacing deltat will be set up in order through 
processing of the above-mentioned step 204,206,207 in the mode "8msec." ->"8msec." ->"8msec." -> "6msec", 
and processing to step 202 - step 212 will be repeated "4 times." and the downstream intake pressure Pm 1 (4) 
in case the counter value i is "4", the downstream intake-air temperatures Tml (4), and those ratios Pml/Tml 
(4) — downstream intake-pressure Pml tfwd after [ this time to ] the prediction period tfwd, downstream intake- 
air temperature Tmltfwd, and a ratio — it will be set up as Pml tfwd/Tmltfwd, respectively. 
[0075] Then, while the counter value i was reset by "0", after discrete spacing deltat is set up equally to the 
above-mentioned predetermined time deltatl (step 216), processing is returned to step 300 shown in drawing 3 . 

[0076] At step 300, the map operation value GA in this time (j) and (it is hereafter written as "GA(0)") are 
computed based on the detecting signal of an air flow meter 41 . Next, based on the map operation value GA in 
this this time (0), and the map operation value GA (j-1) (it is hereafter written as "GA (-1)") computed in the 
last interrupt timing, the throttle-valve through put measured value mt2 in this time (0) is computed using a 
previous formula (24) (step 400). 

[0077] and this throttle- valve through put measured value mt2 (0) — being based — each previous formula (15) - 
(17) — using — the downstream intake pressure Pm after discrete spacing deltat (i), the downstream intake-air 
temperatures Tm (i), and those ratios ~ Pm/Tm (i) is computed, respectively (step 500). In addition, the 
downstream intake pressure Pm computed in this step 500 (i) the downstream intake-air temperatures Tm (i) 
and those ratios — the value computed to this interrupt timing among Pm/Tm (i) with "Pm2 (1)", "Tm2 (1)", and 
"Pm2/Tm2 (1)", respectively It writes "Pm2 (0)", "Tm2 (0)", and "Pm2/Tm2 (0)", respectively. [ the value 
computed to the last interrupt timing ] 

[0078] each following formula (31), (32), and (33) — being based — the downstream intake pressure Pmtfwd 
after [ this time to ] the prediction period tfwd, the downstream intake-air temperatures Tmtfwd, and those 
ratios — while Pmtfwd/Tmtfwd is computed, the inhalation air content mctfwd in a cylinder after the prediction 
period tfwd is computed by being based on the formula (34) obtained from the above-mentioned formula (14) 
(step 600). [ next, ] 
[0079] 

[Equation 23] 
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Pmtfwd = PmMfwd " Pml <•> + Pm2 «» •••(31) 

^_ = ^t±_^± i0) + ^2_ m ... (32) 

'mtfwd I m\tfwd i m\ 1 m2 

^=T^ "' (33) 

F mtjwd 



In this way, after the inhalation air content mctfwd in a cylinder is computed, the initial value of data processing 
in next interrupt timing is set up, respectively (step 700). 

[0080] Drawing 5 is a timing chart which shows the setting-out mode of such initial value. In this drawing, 
while engine operation is started in time of day tl, this routine shall be started. 

[0081] Initial value [ in / as shown in this drawing / the interrupt timing (time of day t2 and t3) of the 2nd 
henceforth ], Namely, the downstream intake pressures [ Pm / Pm (0) and / 2 ] 1 (0), the downstream intake-air 
temperature Tml (0), In Tm(s)2 (0), these ratios Pml/Tml (0), Pm2/Tm2 (0), and a list, the map operation 
value GA (-1) It is set up as follows using the value computed in the last interrupt timing as the values (Pml (1) 
etc.) or the value (GA (0)) in this time after discrete spacing deltat, respectively. 
[0082] 

[Equation 24] 
Pml W «- Pm\V) 

Pm2 <°> <" Pm2 CO 

T m \(o) <- T m \(D 
T m2 i0) 4- T m2 W 

J m2 'ml 
GA{-\) <- GA(Q) 



On the other hand, about the initial value in the first interrupt timing (time of day tl), it is set up as follows from 

the calculation value in the last interrupt timing not existing. 

[0083] 

[Equation 25] 

Pm\ <°> «- Pa 

Pm2 «» <- Pa 

T m \ (0) <- Tha 

Tml (0) *- Tha 

P P 

^(°>-^- 
l m\ i ha 



l ml *ha 
G/4(-l) <r- 0 

Thus, after setting up each initial value in next interrupt timing, processing of this routine is once ended. 
[0084] As explained above, with this operation gestalt, a fluid model is built based on a law of mass 
conservation in the heat energy conservation law list about the inhalation air of the downstream part 28. The 
fluid model, Based on the downstream intake pressure Pm and the downstream intake-air temperature Tm, the 
downstream intake pressure Pmtfsvd after the prediction period tfwd and the downstream intake-air temperature 
Tmtfwd are presumed from this time in the throttle-valve through put mt and the upstream intake-air 
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temperature Tha, and a list. Furthermore, he is trying to predict the inhalation air content mctfwd in a cylinder 
after the prediction period tfwd from this time based on them. 

[0085] (1) Therefore, even if it is the case where the downstream intake-air temperature Tm changes from this 
time transitionally, when the transitional change has been grasped, the inhalation air content mc in a cylinder 
can be predicted, and a very accurate prediction result can be obtained. 

[0086] In case the fluid model of such inhalation air is built, the law of mass conservation and heat energy 
conservation law about inhalation air of the downstream part 28 (2) Moreover, a previous formula (12), He is 
trying to presume the downstream intake pressure Pmtfwd after the prediction period tfwd, and the downstream 
intake-air temperature Tmtfwd from this time by a discrete formula's as shown in (13) expressing, and 
performing processing which repeats the operation based on these dispersion type (12) and (13) the number of 
predetermined times. Therefore, the downstream intake pressure Pmtfwd required for prediction of the 
inhalation air content mctfwd in these cylinders and the downstream intake-air temperature Tmtfwd can be 
presumed now based on a comparatively easy algorithm. 

[0087] (3) In case the above-mentioned repetitive operation is performed a predetermined interrupt period, he is 
trying to compute the downstream intake pressure Pm (0) and the downstream intake-air temperature Tm (0) in 
this time used to next interrupt timing in this interrupt timing furthermore. For this reason, not using the sensor 
for detecting the pressure and temperature in the downstream part 28 etc., ** can also ask for these downstream 
intake pressure Pm (0) and the downstream intake-air temperature Tm (0) serially through the above-mentioned 
repetitive operation, and can attain now simplification of the configuration of inhalation air content prediction 
equipment. 

[0088] (4) By presuming the downstream intake-air temperature Tm in this time (0) especially, the need of 
generally using the low general-purpose temperature sensor of responsibility is lost, and lowering of the 
predictability resulting from the detection error of such a sensor can be avoided now. 

[0089] Moreover, he is trying to calculate the throttle- valve through put estimate mtl presumed from engine 
operational status, such as the throttle opening TA, and the throttle- valve through put measured value mt2 
computed based on the detecting signal of an air flow meter 41 as throttle-valve through put mt with this 
operation gestalt, respectively. 

[0090] Here, when the downstream intake pressures [ Pm / Pm (i), Pm 12 I (0), and / 2 ] 1 (1) called for using 
these throttle- valves through put estimate mtl and the throttle- valve through put measured value mt2 and the 
downstream intake-air temperature Tml (i), Tm2 (0), and Tm2 (1) compare with a actual intake pressure and a 
actual intake-air temperature, it is checked by experiment of this invention persons that there are the following 
inclinations. 

[0091] That is, about the downstream intake pressure Pm 1 (i) called for using the throttle- valve through put 
estimate mtl, the downstream intake-air temperatures Tml (i), and those ratios Pml/Tml (i), there is an 
inclination for the temporal response to reflect the actual temporal response of an intake pressure or an intake- 
air temperature in accuracy especially. On the other hand, about the downstream intake pressures [ Pm / Pm (0) 
and / 2 ] 2 (1) called for using the throttle- valve through put measured value mt2, the downstream intake-air 
temperature Tm2 (0), Tm(s)2 (1) and those ratios Pm2/Tm2 (0), and Pm2/Tm2 (1), there is an inclination for it 
to be well in agreement with an intake pressure, an intake-air temperature, or the actual value of those ratios at 
the time of steady operation. 

[0092] For this reason, if it is at the time of transient operation by computing the downstream intake pressure 
Pmtfwd after the final prediction period tfwd, and the downstream intake-air temperature Tmtfwd using each 
previous formula (31) - (33), while being able to make that calculation reflect the temporal response of a actual 
intake pressure and an intake-air temperature in accuracy, it is at the time of steady operation, [0093]. [ i.e., ] 
[Equation 261 
FmtXjwd = PmlW 



P mt\Jwd = 
T mt\fwd T n 



ml 



m 



It is [0094] when ******(ing). 
[Equation 27] 
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Pmtfwd = Pm2 (°) 
P mtfwd _ P m2 

T mt fwd T m2 



(0) 



Since it becomes, the downstream intake pressure Pmtfwd and the downstream intake-air temperature Tmtfwd 
come to be called for only based on the throttle- valve through put measured value mt2. 
[0095] (5) therefore — according to this operation gestalt — the time of steady operation — be — the time of 
transient operation — be — the downstream intake pressure Pmtfwd after the prediction period tfwd and the 
downstream intake-air temperature Tmtfwd can be presumed from this time to accuracy, as a result the 
inhalation air content mctfWd in a cylinder can be predicted more now to accuracy. 

[0096] (6) Furthermore, with this operation gestalt, since the thermal equilibrium model about the glass layer 
414 of an air flow meter 41 is built and he is trying to calculate the above-mentioned throttle-valve through put 
measured value mt2 based on the reverse model, the response delay about the temperature change of this layer 
414 can be foreseen beforehand, and it can ask for this. Therefore, the more exact throttle-valve through put 
measured value mt2 can be calculated, and the downstream intake pressure Pmtfwd and the downstream intake- 
air temperature Tmtfwd can be further presumed now to accuracy based on this. 

[0097] The 2nd operation gestalt of this invention is explained focusing on a point of difference with the 1 st 
operation gestalt below [the 2nd operation gestalt]. 

[0098] With this operation gestalt, the point he is trying to apply the inhalation air content prediction equipment 
concerning this invention to the internal combustion engine which has an intake-pressure sensor as a sensor 
which detects an inhalation air content is different from the 1st operation gestalt. As shown in drawing 6 , this 
intake-pressure sensor 46 is attached in the surge tank 22, and detects the pressure Pm of the inhalation air of 
that interior, i.e., a downstream intake pressure. 

[0099] Next, the detail of the prediction processing by the equipment of this operation gestalt is explained. 
Drawing 7 is a flow chart which shows an example of the procedure at the time of predicting the inhalation air 
content mc in a cylinder, this processing of a series of is performed by the arithmetic unit 30 the interrupt period 
(for example, — "— every 8msec") for every predetermined time. In addition, at step 100 shown in this drawing 
7 , and step 200, since the same processing as step 100 and steps 200-216 which are shown in previous drawin g 
3 and previous drawing 4 is performed, explanation of the content of processing is omitted. 
[0100] After each processing of these steps 100,200 is performed, at step 350, the downstream intake pressure 
Pm in this time (i) and (it is hereafter written as lf Pm3(0) M ) are computed based on the detecting signal of the 
intake-pressure sensor 46. Next, the downstream intake pressure Pm detected in the downstream intake pressure 
Pm 3 (0) detected in this interrupt timing, and the last interrupt timing (i-1) (it is hereafter written as "Pm3 (- 
1)") Based on the degree type (35) obtained from a previous formula (15), the throttle-valve through put (to 
namely, the last interrupt timing) mt (i-1) (henceforth "the throttle-valve through put measured value mt3 (-1)") 
in front of discrete spacing deltat is computed from this time (step 450). 
[0101] 

[Equation 28] 

m/3( .,> = _i_j ^-JW-.> .JL_ + To{a . Pm m _,)}.. . (55) 

Further And this throttle-valve through put measured value mt3 (-1) and a previous formula (16), It is based on 
the following formula (36) obtained from (17), and (37). The downstream intake-air temperature Tm in this 
time (0) The ratios Pm3/Tm3 (0) of the downstream intake pressure Pm 3 (0) and the downstream intake-air 
temperature Tm3 (0) are computed by the list ("Tm3(0) [ list ]" is written hereafter), respectively (step 550). 
[0102] 

[Equation 29] 

-Jsi (0) = <-D + At . £ - [m t2 (-i) - — &- - (a • /^ 3 <-i) - bj\ ■ • -(36) 

T m3 T m3 V \ *m3<-» ' 

r«3<o> = -£^ -(37) 

http://www4.ipdl.ncipi.go.jp/cgi-bin/tran_web_cgi_ejje 1 0/1 7/2005 



JP,2001-041095,A [DETAILED DESCRIPTION] 



Page 14 of 15 



In addition, in the above-mentioned step 450 and step 550, when computing the downstream intake-air 
temperature Tm3 (0) etc. based on a formula (35) and (36), the value which the value into which the upstream 
intake-air temperature Tha was read to the last interrupt timing was used, and was set up based on the last 
engine engine speed NE and the valve timing VT in interrupt timing about each constants a and b is used. 
Moreover, when this time is the first interrupt timing, the value in this time [ intake-air temperature / Tha / 
upstream ] and the value to which each constants a and b are set based on the engine engine speed NE (= "0") 
and valve timing VT at the time of an engine halt are used. 

[0103] Next, the formula which permuted "Pm2 (0)" and "Pm2/Tm2 (0)" by "Pm3 (0)" and "Pm3/Tm3 (0)" in 
each previous type (31) and (32), respectively, each formula (33) and (34) — using — the downstream intake 
pressure Pmtfwd after [ this time to ] the prediction period tfwd, the downstream intake-air temperatures 
Tmtfwd, and those ratios, while Pmtfwd/Tmtfwd is computed The inhalation air content mctfwd in a cylinder 
after the prediction period tfwd is computed (step 650). 

[0104] In this way, after the inhalation air content mctfwd in a cylinder is computed, the initial value of data 
processing in next interrupt timing is set up, respectively (step 750). Drawing 8 is a timing chart which shows 
the setting-out mode of such initial value. In this drawing, while engine operation is started in time of day tl , 
this routine shall be started. 

[0105] Initial value [ in / as shown in this drawing / the interrupt timing (time of day t2 and t3) of the 2nd 
henceforth ], In the downstream intake pressures [ Pm / Pm (0) and / 3 (-1) ] 1, the downstream intake-air 
temperature Tml (0), Tm3 (-1), and a list, namely, these ratios Pml/Tml (0) and Pm3/Tm3 (-1) It is set up as 
follows using the values (Pm3 (0), Tm3 (0), etc.) calculated in the last interrupt timing as the values (Pml (1) 
etc.) or the value in this time after discrete spacing deltat, respectively. 
[0106] 

[Equation 30] 

Pm\ <0) <" Pm\ (0 
Pi«3(-0 <- Pm3m 
T m \ (0) - T m \ (I) 

J m\ 'ml 

1 m3 J m3 
Pml(-» <" Pm3 (°) 

On the other hand, about the initial value in the first interrupt timing (time of day tl), it is set up as follows from 

the last calculation value or last detection value in interrupt timing not existing. 

[0107] 

[Equation 31] 
Pm\ (o) «- P a 

Pml ("I) Pa 
T m \ <°> <- Tha 

P m \ P a 
2 m\ i ha 



ml l ha 
Pml t-0 <- Pa 

Thus, after setting up each initial value in next interrupt timing, processing of this routine is once ended. 
[0108] Also in the 2nd operation gestalt explained above, the operation effectiveness of an abbreviation EQC 
can be acquired about (1) indicated in the 1st operation gestalt, (2), and (4). Moreover, about the operation 
effectiveness of (3) indicated in the 1st operation gestalt, since the sensor which an air flow meter 41 becomes 
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unnecessary at the reverse of that for which the intake-pressure sensor 46 for detecting the downstream intake 
pressure Pm is needed separately, and detects the downstream intake-air temperature Tm at least is not needed, 
at this point, simplification of the configuration of inhalation air content prediction equipment can be attained. 
[0109] Furthermore, the downstream intake pressures [ Pm / Pm (0) and / 3 ] 3 (1) called for based on the 
detecting signal of the intake-pressure sensor 46, the downstream intake-air temperature Tm3 (0), Tm(s)3 (1) 
and those ratios Pm3/Tm3 (0), and Pm3/Tm3 (1) have the inclination for it to be well in agreement with an 
intake pressure, an intake-air temperature, or the actual value of those ratios at the time of steady operation. 
Therefore, the operation effectiveness equivalent to (5) indicated in the 1 st operation gestalt also according to 
this operation gestalt can be done so. 

[0110] Each operation gestalt explained above can also change and carry out a configuration as follows. 

- With each above-mentioned operation gestalt, as an amount of the inhalation air which passes a throttle valve 

24, although the estimate based on engine operational status, such as the throttle opening TA, and the actual 

measurement based on the detecting signal of an air flow meter 41 or the sensor of intake-pressure sensor 46 

grade are calculated, respectively and the inhalation air content mc in a cylinder was predicted using each 

[ these ] value, based on either these estimate or an actual measurement, it may be made to perform this 

prediction. 

[01 1 1] - Although it was made to predict the inhalation air content mc in a cylinder with each above-mentioned 
operation gestalt supposing the internal combustion engine 1 0 having the device in which the valve timing of an 
inlet valve 26 is changed, also in the internal combustion engine which does not have this device, the inhalation 
air content mc in a cylinder can be predicted according to the same procedure. 

[01 12] - When the device which controls the flow of the inhalation air of the inhalation-of-air path 20, for 
example, the swirl valve which controls the reinforcement of the swirl generated in a combustion chamber 29, 
performs the above-mentioned prediction again in the internal combustion engine formed in this inhalation-of- 
air path 20, set up the above-mentioned constants a and b, the downstream intake pressure PmTA at the time of 
steady operation, and the throttle-valve through put mtTA as a function of the controlled variables (for example, 
opening of a swirl valve etc.) of the above-mentioned device. 

[01 13] - In prediction processing of each above-mentioned operation gestalt, although the empirical formula 
(14) which makes this air content mc the function of the downstream intake-air temperature Tm and the 
downstream intake pressure Pm was used when calculating the inhalation air content mc in a cylinder, the 
method of calculating this air content mc is not limited to this. 

[01 14] - Although the downstream intake pressure Pm and the downstream intake-air temperature Tm were 
computed with each above-mentioned operation gestalt by discretizing each above-mentioned formula (8) 
obtained as a fluid model of inhalation air, and (9), and performing a loop operation, you may make it compute 
these downstream intake pressure Pm and the downstream intake-air temperature Tm based on the operation 
expression which solved each [ these ] type (8) and (9) analytically. 

[01 15] - Although the above-mentioned prediction period tfwd was set as the value suitable for calculation of 
this fuel oil consumption with each above-mentioned operation gestalt that the inhalation air content mc in a 
cylinder used for calculation of fuel oil consumption should be predicted, this prediction processing is 
applicable also to other engine control, such as ignition-timing control, by setting up the prediction period tfwd 
suitably. 

[01 16] - Although the interrupt period of a prediction manipulation routine shown in drawing 3 or drawing 7 
and discrete spacing deltat were equally set up with each above-mentioned operation gestalt, this discrete 
spacing deltat can also be changed "1/2 of range which can determine initial value in next interrupt timing, i.e., 
interrupt period," twice etc., for example so that it may become 1/n time (n: integer). 
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[0 04 1] 
[f 9) 

±3t (14) lC*SW<\ Taj % Tbj ttV^nt^R 
0') = P-tf-D + A/-*- * 



tt^^HS^^J:«9S^5>tL. SKJHSK3 o<0**y 3 2 
IcHft^s/ri LTl£t6*tLT^5o r Toj 

10 KAffi«OSfllJftK "Cfc*. 

[00 4 2] (1 4) .&m^*ks .hiE&it 

(12). (13) It (15), (16) (0 

[0 0 4 3] 



[»1 0] 



F ■(«,,<,-!) 

P P R ( 

jr- CO = y 2 - 0 ~ i) + A/ ■ — • j^m,(i - D - — 



/jo ~ To 
To 



■ t) 



•(15) 



TS6ffl«SUBTmfl* ZtLb^ (15), (1 
6) a*&»P>*x5«T<&5£ (1 7) Srffl^T^SZi: 

[0 0 4 4] 
[mill 



-••(17) 



PJ 



:m (is) - (17) ^"e>w?>^<cJ: a 

ftfiTO^nyh/U#iiSlmt, T86(R!ia«U3EPm, T 
SEffllffiSliBTnu 25it5±«EftflK«?aTha»cS<5^^T, §1 

OT^KaiEPm^T«iRflK^TraSr«l^i-5wi: 

T8ltflBfa«faTinS:S: (14) icttAtSwilcJ:!}, 

[0045] ^/c. *j8jffijga8<o^»j«ia-ctt, ±e* 

[0 0 4 6] 4i\ x77o-^^4 l*>»tHflr»fc 

4 i 2^&»fcixSl!fta^KA2«4lcjS;i:rS{bi-5 
-^-^41 ottffl«*t KA35*fi t <OB8{»tt^«)* 



30 



40 



[0 0 4 7] La>U>WS?> % JilS^7^^4 1 4(1, ik 
j^4 l 2 ttttfcUT— ttfc:SaSftaUc&v*fc«>* ^ora 
J2Tgsf±«a5ffl4 1 0Oi£«*iBiS+S»A£«tf>»<0 

Affift »a*SgflS L v ^* SlteW i:J)oTll JilSSS 
y ^ fe*«> &*u5BlAffl**ttH<OSRAffiK*i: - 
SrTSfcoo. HKAffi«ft^ffi-ffc-*-5iS«ai«5B*Jwfo 
oTtt, H»^s/^»b*«>e>ixSKAffi««:^Jt^K 

[0 0 4 8] *ZX\ **»»tt-Ctt, ^ ? Lfc#7* 
Mai 4^^^b«coVN-c^>Jc£:^51ixSr J f^^aUT 
I^SA^M**^^^!:, |5]J§4 l 4dfg-r3£* 
^i^T/^tHL, *0>i^f*/MCg<5^TKAffi* 

0»AS«tS:^D 5/ h/U#ilias:mtt U x77n 

So 

[0 0 4 9] % t&»4 1 2<Or£./££ TThtJ 
t (i2#I) % 1 2a>fe#7*JB4 1 4IC&X 

SitSSftfti:* fU^4 1 4^bSRA3S«lceitSixa» 
iiOft^l^Pj, X=7*m 4 1 4CD?&STgstf>^fc 
id (Tgs) /d trtt, KTOS (i 8) oJ:?i:»t 

[0 0 5 0] 
[»1 2] 



(7) 
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12 



* 27 T gs = b • (77,, - 7>) - (c + . V^) . ( 77,, - r*a) ■ ■ -0 «) 



taj , tbj . rcj t tdj fft»4 i 

t«ij»4 1 2 <tomO«ReJt^ ^7^g4 14t8A 

[0 0 5 1] £1RiMBR*(Ctt, ^/7^I4 1 4 

Sfc£>, ftS(19), (2 0) {o^fj: 5t (l 
8) ^&522, EP*>^7^^4 1 40i&ff^{fcfid (Tg 
s) /d tl± fOj \Z*9 % ▼y^SWGAt 
*n h/U#ifiifl4mti:fi«U<ft5. 

[0 0 5 2] 

Itl3] 

o = ^-(r A ,-r g ,)-(c+o-^-(7-M-rA fl ) "<h>) 

GA= mf ---(20) 



i) wens, 

[0 0 5 3] 
[mi 4] 



(19). (2 0) frhVXO*. (2 



••(21) 



HI-, rco^; (2 1) Sr/flv*T5t (18) Srg^-T* 
tf^xm 4 1 4<DiSfiTgsJ&s?H*£*LfcftS; (2 

2) aq#<=>n5. 

[0 0 5 4 J 
[ftl5] 



•(22) 



1^ 



10 



20 



30 



•(23) 



rost (2 2) oaEffltt. BStkBaPHSr t j 
-rsir&xt (2 3) <D£o\z&-fzktfX'ZZ> 0 
[0 0 5 5] 
ffil 6] 

dyfCA _ yjGAijy - jG/fr/ - I) 
~ A/ 

feoT, ±K** (2 2) . (2 3) *»fc«T«>Blft5e 
(2 4) 3&S«*aSixS. (2 4) fclteV* 

•Sct'fc 9 , *fc, tsarto r j j tt^Eifi, r j — i j 
ttOT IsHK £ * *t LTV>6e 

[0 0 5 6] 

[ftl 7] 

Wf UI-(V^^^ ^"^- ,> 4 1 -(24) 



rcoiS (2 4) *»&9it>ft»ttj:5f;: % TOt'O^ny 
(j)*. MttHHIA tmo-^y^*|MtGA(j-i)i:lc:B 

[0 0 5 7] y h/UBajffiTA«fO«B8iBett 

(D&jK. (2 5), (2 6) lCj:c.Tfl$n5o 
[00 5 8] 

mis) 



>*{P m ,P a ) "-(25) 



— If— 1 



2/r 



2 



a a U+J 



AT— I 



••(26) 



wwT% r^j litl)7^7 (*n h/u# 2 4<£>(J3 a [0 0 6 0] 

ogg*«&. rpaj fi*mm (-jmo -c*> [fti 9] 

5 0 -£ti. TAtj (i, ^D3/h;u#2 40BDi«-e 
ZtMxu y h/^Ba^TA^e>-SS6<j{c*^-5w 

[0 0 5 9] ^n^/uKTA^-SCi^ 
Sft3S*iBtell*«rfflS1-5i:, 5C (2 5) ttfcsC (2 

7) ^ijicfts. .so 



ha 



^r^o) ^(27) 



fflmtTA2a.TJ { Ti?pLfll"IK«/EPmTAfi, JiSSffiilffltf* 2 7rt>*b 



(8) 
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*<TA = T 3 -' ( a * F mTA ~ *) •■•(29) 
* ha 



*d y h/^iffiiSfimtTA^^(iTffi{R»JK^i±PniTA<t. 

[0 0 6 1] (27) Si (25) li 10 

(2 8) (DX :[:ttlt^T^6o 
[0 0 6 2] 
[t20] 

Ai1Brt«ASft*mc£aS»L< 4 9 (mtTA = mc) , 
T»«yK»ffl.Tmi:JtatiW«««aThai:dS*L< 
KZtzfr (Tm=Tha) , ±155*; (14) ttft^ (2 20 

9) <DJ: b&T'ZZ. 

[0 0 6 3] 
[§£2 1 ] 

^LT, rcO^: (2 9) *m^Zh. JzlSxt (2 8) 

JEKUfcst (3 0) (OX 5lw*tri#-ct«. 

[00 6 4] 
182 2] 

-•e •(-".»-») •••<*» 

(30) ^e>n&^«j: 51-, 

JEPnu *:^J£Pa, 3fetffC5^jIiE8#T*fc5 £ fi£tfc 
<t #<OT^RiJB&^/EPm<^MP ra TA (TA, NE, VT) 

[0 0 6 5] ^Hlg^lT'ii, :^J: 9twS-^Jo^fe»c 

fc&zfc (14) - (1 7) <h{CS^^T!EB$^^^^ 40 

mt^Q y K/U#igiga8!!j^«i:mt2i LTES'if Co 
[0 0 6 6] *|t*fiJgaS^|$EtcJ:-5^»l«i3ffio 
Baotcou^TEI3 3lVH4^fif-ti:#IHLTttMi-5. H 
3Wi4li, «rt®A25«fflmcS:^8IIJ1--5l»^^)lR 



T 8 msec. J g) "CH?r£ix5o 

[0 0 6 7] io»ll:BLTrt, *i\ 8ftAffi**«r 

o) o zco^mmmtM^ m^com^mxcomnm 

£ fiS'J ojRjpfpAtt S «ka «: ii IS r K Ui $ tt 5 1 
£rii lSTKUS*ix5ft»r^«rtaAS»ftmclca-^v^ 

[0068] jsfc(c m$&frb^mMmfwd&<DTWiM 

»«UEPnu TSMRdKflHaTnu &U^;h,<bOtt: Pm/Tm 

it* (77^/2 0 0) 0 r*h Jt^TT'fi. ^n^h/U# 

OT«WJ8R«JEPm, T«E«K«aTm. Xtf-tntotf: 
Pm/TmSr-tfL-piV rpmltfwdj „ TTmUfwdJ . 

rPmltfwd/Tmltfwdj k % •£0#fcBiS + -COT 

ScMKSaJEPnu T«{M»«&Tiiu Sr^-tn^^JtPm 
/TmS:, rp ra l(i)j % rTml(i)j % TPml/Tral 
(i)j ir^rix^tu^iSi-So 

[0 0 6 9] wOT8IEffilK«ffiPinltfwd&t/T8E«J®« 
?aTmltfwdO»tfc5{C^LTtt, ±fS#5t (15) — (1 
7) Xtf* (3 0) lw»<5<«l|[fc3fSEI»«lt»9 3S 

[0 0 7 0] EP^>. IB 4 5 |C % £-f. 

>^ y > > h£*x (^77^2 0 2) , ftK^SMWBfltf 

>owa«ast«u< (M;tr£ rsmsec.j ) 

tl^o -t LT, ^»J«Brp1tfwd*SwOBf^B#MAtlEJL± 
tfcSifllUfS^S^ (^X^^2 0 4 : YE S) , 
JtE#a (15). (16) OBKMIHA t**J!fJB»IB 
A tl<t^L< I3:^$^L6 (^75/^2 0 6) 0 
[0 0 7 1 ] ^ (30) Clo^t^D v 

(15) - (17) ^cS^v^T, TSKiWKftlEPnil 

(i), TSE-ftiiK«faTnii(i). &v*ne>oit p mi (0/ 

Tml(i)3dS**XrPixKaS*x5 (^fy^2 0 8) o ^ 

it, afto^»]j«mtfwd^e>Bf*B*raAti^«i(s 

tt, -ttf>«JSffi (tfwd-Atl) ^rrfc^^ifflJWrBltfwdi: 
LTlx^t^no (^^^7*2 1 0) o 

[0 0 7 2] z <omizte^®m ffltfwdtf roj » 

JBMtfwd** roj i: 9fi^£4Wr£ftfc*J^ (^^^/^* 
2 1 2 : NO) iCii. 7f^2 0 2 JWReoj&3//pffl[ 



(9) 



15 



tfwd*saf«B*IB A tl*jBT*>4 i: *J»rS;h,5 fc (^fy 
/2 0 4: NO) , RtlKmHSA t *s?jHJO)Rntfwd « L 
<IS:S^iXT (7f^2 0 7) % *7-*y72QS$m 

[0 0 7 3] -t IT, 7^77^2 0 2-2 1 2Wftl3i J 

«k9ES;h,SwifcJ;9. ^mnmtfwd^ roj J^Ttc 

^^><t (^7y7'2 1 2 : YES) % TffiWJiSUE Pml 
(ih Tiffi««»«ffiTnl(ih XtfJt^&OPml/Tml 

^»JKBIBtfwda<OT«lE{MKft£EPinltfwd, TiiS 10 
{Ri]9£gtfSTmltfwd. &tntPmltfwd/Tmltfwd£ Ltt 

;h,^jh4£J£$*L5 U7^2 14) D 

[0 0 7 4] 0J*.tf. ^^ffli&JflJtflj&S TSmse 

c.j Xth*) . 9t<0*TvZfl O.OO&SSra&T^ajJffl 
MtfwdtfS r 3 0msec. j lCgfcj£ § fr*:^, REftMlffiA 

1^112^7^/2 0 4, 2 0 6, 2 0 7^15:iii: 
T\ r 8 msec. J — > r 8 msec. J — » r 8 msec. J — ► T 6 rase 
c.j i:V^o^:®.«T^(cS:^$tt. ^7^2 0 2^ 

KZo ^tt, r 4 j T*^5i:tcOTc?t 20 

flUJKSUEPml (4) , T«HM»«?fiTiiil (4) v Rt*fc 
tb^Oit Pml/Tml (4) an#jfe^P,^j|BjH|Btfwd 
^OT^«^ffiPmltfwd. T«{R(Sft»ffiTmltfwd* & 
tfkt P ml tf wd/ Traltfwd t l^X^t^tltx^ ZtlZ^k 

[0075] ^ot, ^^v^fir i ^ roj icy-tr^ h 

IXKSSixfcfS (77^2 l 6) % MIHH3{:^ 
t^ry/3 0 0 (cjp|£*L3 0 
[0 0 7 6] ^77/3 0 0Tli, m^^a— > — ^4 30 

(j) (J^T, TGA(0)j t»«Bi"5) a*lfai3*L*. ft 
(c, ^ co^^t'(Dv . 7 G A (0) <!: ^0^®^ 

^ s >^ic:fc^-cS¥ai $ it*:^ * ^»*ttG a (j-i) 
rGA(-l)j tm?rfZ>) t\zm<$$, 9t<r>H 

(2 4) &m\,^Tm$&X'<D*v v h/^iliSfigB^e 
mt2(0)^|fUiSnS ^7^4 0 0) . 

[0 0 7 7] ^U, ZO-Xn y h/U#iiiSftai^fflm 
t2(0)tc£*3#. ^cO#5t (15) - (17) Srffll^T 
BtfiJcP»1fSA tm'OTSSmi|55«EEPm(i), TiStt«5KSW&T 40 
m(iK ^.I/ j ett^COM:Pm/Tm(i)^^tL-eiX^ttS^ix 
5 (^7^5 0 0) e ft. I^ry/SOOl:^^ 
T^£;h,£T^WJ@mi£Pm(ih TfiKflBKSUaTin 
(iK XU<*;h,&C0ttPm/Tm(i)O5*>. ^[EjOtiJiA* 

-< ^v^TSttJSixSfflSr-tix^tt rp m 2(i) j , rim 

2(1) J . rp ra 2/Tm2(l)j £, mJlnKafiJiA* -< S >^ 
T'SLtB^tL^lfilSr-tn^n rp ; ..2(0)j , TTm2(0)j % 
. rp m 2/Tm2(0)j <h^f£-f3 e 
[0 0 7 8] J^TO#xC (3 1). (32). 

(33) \cm^^xmn^h^Biiw^fw6m(OT^m so 
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KSOiPmtfwd. Tffi«4!a«rfiTintfwd* RV^tlbOtt 
Pmtfvtd/TmtfwdtmttiZlriZtthlZ^ ±K& (1 
4) frb'ftbtlZX. (3 4) lwS^T^»JIWMtfwd« 
tc*5W5lSjrt»A3£ftftmctfwd*sjfUiSix5 (*7* 
76 0 0) 0 

[0 0 7 9] 

[f23] 

Pmtfwd = ^mlf/v^ - P«l (°> + P*>2 «» 
P mtfwd _ P m\tf«d P m \ [ P m2 



l mtfwd 



TmltfWd 
_ ?mifwd 
P m//wd 



••(33) 



(0) 



••(31) 

•••(32) 



'mtfwd 
P mtfwd 



n ctfwd " 



P mtfwd 



r p UTlBrt8RAffi«*mctfwd36SWB**L^:«. .&[§] 
(^7y/7 0 0) o 

[0080] D5ti. ^tzftmmcD&Tzmm&^-t 

[0 0 8 1 ] Hiatc^-rJ: ?IC. 2 Eh] @ J^^fiHi* >f 
(B##Jt 2. t 3) (C*5*t««»H«. gp*>TSEiffl 
P^^/±Pml(0). Pm2(0). TiSft^JeSMtaTml (0) . Tm2 
(0). r^^(OltPml/Tml(0). Pm2/Tm2(0), ^t>* 
^y/ffltJHKGA(-l)(i, HUlE]OW5i^>f Si-im** 
V^T^Pp E lftl5A tM (Pml (1) m m 
TOI (GA(0)) t LTffffiSn^filSrfflV^T-tix-e 

[0 0 8 2] 
115:2 4] 

/'ml (0) Pml (0 
Pm2 <°> *" Pm2 0) 

7'mt(0) «- r ro i(D 
r m2 1°) T m2 CO 

(") 



(1) 





(0) «- 










(0) «- 




7"m2 


r„ 2 



[0 0 8 3] 
[£C2 5} 



(10) 



17 

Pm\ W - Pa 
Pml (°) - Pa 
T„\ <0) - Tha 
Tml (0) - Tha 
P„ 



'ml J ha 



A. 

[0 0 8 4] *jttS«ffi-Ctt. T 

^r/^t, h/u#iiigfl:mtXV±a£«R«aTh 
a v ttW^T««<»JKffiPiiiftt/T««l»*aTnii: leg 
^V^tWAi»& ^ffl'JJW BQ t f wd'fg <D TiJfciR'J aa«UE P mtf w 
dXu«TiK««ft«aTintf wd«r*6£ U 3E(e-tn 6 

d£^$0i-5 «fc ?{cLT^«5 0 20 
[0 0 8 5] (1) Sot, TXMKXfflTmtfattjft 

[0 0 8 6] (2) :5lfc8A$M**7 

KftftffflURt/Sl^i^^^ftffliJSr^O^; (12), 
(13) »c^Six5J:5fc«K*i;Lj:oT*au cix 
e>REi&a (12) . (13) lca^<WW«5rJ£E!Ri* 30 

tfwdSOTi«W»«JEPmtfwd&t/TaElR!l««aTmtfwd 
^t5<t9(CLTV^ 0 ^ot, :i;h,<bfiJl*lRA£ 
^ftmctfwdo^ffl*c^B<cT8E««K»ffi PmtfwdXt/T 
8K(!l!lRftaTintfwdS:ltttMf»**r/^^!; XACI-^S 

[0 0 8 7] (3) JEIC % ±EiHkigteaS:Bf*OWi&ffl 

an** r-co tsk»]« sue p m (o) a T8E«iKsea t m (o) 

o 0 Z<Dtztf>. TSc<M*!5#2 8 rtOEA-^lKft&ttm-*- 

m(0)XVT^«K«f&Tn(0)fr±E|ft3S*aaS:aDDTS 

flsSrEl 5 C £ *«-e * 5 J: 5 
[0 0 8 8] (4) aB#j£TWT««»*»Tin 
(0)«r*£1-*::i:fc:J:'3. -«»wJC*ttOffiV^iflfflO 
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[0089] *sa6»«-ett, ^o^h/u^aia 

h/^#iligaa'l^flSmt 2 £ £tixt';ft#tf>S cfc 9 J- L 

[0 0 9 0] ~ ~T\ ZfLbX* v h/U#iiiS*t£Sffl 
mtlR^ny h^#iiia£SB'l5£Mmt2S:fflV^T3ft>!) 
fe*lST8ttft<»«cUEPnil(i), Pm2(0), Pm2(l)&tfT 
««»MfiTml(i). Tm2(0) f Tn2(l) tt % XRO«X 

[0 0 9 1 ] BP*> % ^cy h/U#iiifi*«|£fimtlS:ffl 
V*T#a& fett5 T8MW5RSUE P ml (i) . TaMBSftStf&Tm 
Hi). .St^-tn&OitPml/TmKOlCiaLTtt. » 

fflSSJSffimt 2 SrfflV^T**P>ix5TiK«JK«iEPm2 
(0). Pm2(lK -TflEftl|KKaTin2(0). Tm2(l)« X.tf-t 
tue>COJtPm2/Tm2(0). Pm2/Tm2 (l) KPJ UTI*. 

*aME*»c*jir ^TKSUE-^RStffi* Hfcv ^rt:t*^ <b olfco 
JtB©«i J: < -iW--6«ia]*sfcs 0 

[0 0 9 2] rCD7t#>, (3 1) ~ (3 3) 2r 

ffiv^r*»w*^M»PBitfwdaoTa£{a(»«tjaEPnitfwd 
ar/TSsaflKSCfaTmtfwdSriwi-r 5 t ic * 9 % 

fc«H»lc-tO*UiKRI*S-fra - £ *s-e# S 4: £ fc 

[0 0 9 3 ] 
[»2 6] 

PmtXfwd = PmL(°) 
P mt\Jwd _ P m \ 



(0) 



l mtljwd 



[0 0 9 4] 
[ft 2 7] 

Pmf>rf = Pm2 CO) 
P mifwd _ P m 2 
T mt fwd 



(0) 



iWlt^fe, T*«BR»E Pmtfwditf T««l»* 
a T m t f wd^S ^ n y h ifi Sffl itfflC m t 2 (DfyltZ. 

[0 0 9 5] (5) t^oT. *S«fi»aWCJ:*Ltf. 

wdaOT^fflllR«JEPmtfwdXVT^B«a«iaTmtfwdS: 

[0 0 9 6 ] (6) 35{C X #H*£JI2!gt*te, ^77o- 



(11) 
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^4 1^7^14 i 4\^m-t^^w-m^r : f/^m 

«Jj£ffimt2 3r#£>3 X 0 \ZL,X\,^Z£lZ> s [1114 14 

h/u#iiififl:»JSffl[mt2fe*«>5wi:^-et. wft 
KS^TTSEilUIKSUE Pmtfwd^TattflKftJ&Tmtfwd 
^«^€:5l-iE5tlic^ 5 ^ t 5 <fc o left 5. 
[0 0 9 7] [ff 2 <DM1&Bm £*T\ *3g9iOSB 2 O 

Hi-fi. 

[00 98] *Ha£»J6-Ctt. «A£ftftS:#ttSi-3-fc 

fiRAffifci^M&BfcBffli-* <fc 5 LTir^* jft*sjB 
1 OHjKJBffi^ttiSLTV^. w608S^EEir^i^4 6 
(A. {H6(c^iTJ: 5 -j— v?* v-^ 2 2(cRt? ft»+b 

JSEPm€r«ffii--5. 

[0099] *mw&m<ommz&z>^m%ii&<r> 

(-D 



20 



i&JSJW (fflaitf rsmsec. j ft) T*H?f$ix5o ft, 13 
[H ^"1 0 0S^7y^2 0 OT'ft, $fc 

O0 3Xt/El4 tw^-J-^T-y^l 0 0^T77'2 0 

o — 2 1 6 km&^mmt-tmzftzizito. tcD&mft 

[0100] Cti^^f^/lOO, 2 0 0O4ftf^ 
• ^T^JtfcS, ^^5/^3 5 OT'fi. ffi«JE-fe>*4 6 

^*tH«*fc*^l^T»B»jaTOT««R«EPni(i) 
io (j^T, rp m 3(o)j t*Ef5) ^WtUSixS. jfc 

WKSC£EPm3(0)^t^SHs]OS»Iii^-f ^ ^iCfcJ^T^ 
UiSixfcTiKffliJKStEEPniti-i) (»T. rp m 3(-i)j t 
*E+S) ^fecoxt (15) ^feWb^Sft* (3 

5) irjcso-v>r. %nFja^&Ri«mHKA tmr-co (hp 

(i-1) (^T. r^n y h/u#ffiigM]5£ffimt3<-l)j 

[0101] 
20 [§C2 8] 



7*» 1 At 



A/ 

^tt, JEfclwO^n y h^#iiigapJ^Mmt3(-l) 

5feO^ (1 6) , (1 7) a>^#P>;ft,5&0:i£ (3 
6) , (3 7) klzm<5\t^T. W%&X~<DTffi®mf5.& 
Tm(0) (J£*T\ rrra3(0)j <t^i5-r 5) . El^TSE * 



JTTyf + To ( a 'PmW -*)} "-(35) 



■J 52 - CO) = ia- (-D + A/ • A • L, 3 (-., - — 3> ( 



7 m3 



^«IK*ffiPm3(0) tTJ«mfl!»ftffiTm3(0) t^Jt Pm3/T 
m3(0)^^tLjtUJ*jXS Uf^5 5 0) 0 
[0102] 
[*29] 

^.(.., ...(36) 
(37) 



ft. ±f£^'7"^^ , 4 5 O^^f yZfS S.OlCjfcl^T. 
« (3 5) Sit/ (3 6) ^S^^TTSE«IKfi«g.Tm3 

(w^srjtaji-sa^ ±8E{Rij»«iaThai*afna^W5& 
a, btcou^rttw®oW5&^>r ^^^T*o«Ba[5]te» 

a. b (i«H»jfcB#o«B3ig4E»N e (= roj ) Xtf 

[0 10 3] ftJC, 9t<D&3Z (3 1), (3 2) twio^ 
rp m 2(0)j % rp m 2/Tm2(0)j rp m 3 
(0)j . rp m 3/Tm3(0)j KEftLfc^i* ^g-xC (3 
3) . (3 4) fc&fflVT, 3RD^ft^^^ffliJ»I|atfwdS 
<OTSEfiHK«/EPmtfwd, T8E<WSR«?aTmtfwd. 
^«bc7)i:kPmtfwd/Tmtfwd^|?£ti^^oi: <b ttw, ^fgiij 50 



JHMtfwdStcfc»t5f8rtKAffi«ftmctfwd^»aiStt 
5 (^7^^^6 5 0) o 

[0104] ^ 5 LTWrt«A£«ftmctfwdas»ttJ£ 

[0 1 0 5 ] leJEllc^l-J: 5ic % 2[H]gJ[^|^coSijiA^^ 
(B*»Jt 2, t 3) lw*5»t5|0J«lit, iP*>TSEflW 
KJWEPml(O). Pra3(-1) % T««»««.Tml (0) . Tm3 
(-1). &Zf\CZtlb<DttPml/Tml(Q) t Pnj3/Tm3(- 

fs (Pmi (i) m i&\,v*%im*X'<omtLx#&bn 

fcffl (Pm3(0). Tm3(0)^) Srffl ^T^iimWTO <t 
[0106] 
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[S3 1] 

Pm\ C°> «- />„ 
^3 (-') «- Pa 

^3 C-l) - Tha 

^ll ( o) «_ A. . 
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Pml C-0 4- 

[0 10 8] «±tt^L3t»2^3llf6^ffi^fcV^Tfc % 
IB l o*i£^«8»ciI5«6 Ufc (l) , (2) % St/ (4) 

m 1 OHJgJgffilClSlffi Lfc (3) OflUBSWfcttloi^T 
14, Tail0Jfa^Pn*«Ui-r5^*O5K«LJEir>'l^4 
6 4i8Ujfcg.Sl;:tt& tOO, i^fc^T^D— ^4 1 
tt^WcafcD; '>4<itT«Bll*«lTi*ft 

[0 10 9] SgEt>t4 6 C0^£Uff-§-|Cg-3 

^T**5>ix5T8£{a!l^/EPm3(0), Pm3 (1) , T 
SKMK«faTm3(0). Tm3 (1) , SU^rtl, OJt P m3 
/Tm3(0) f Pm3/Tm3 (l) 14, ^«51(EB#(i*JV^T 
KjKJE-^SflS« ^v^4^:ixf>0)tO|IBgofli[i!:^< — 

coj«S»«Jci3V^TE«Lfc (5) im^^M^ 

[0 1 10] «±KWLfc*Hte7BiSf4, WToj:9{c: 

• ±E#IIJ!i£7KfigT-(4, ^ d y h/u# 2 4 £2Biai"5» 



>-^4 6 ^coir>iho^^{f-^Jc:So'< ^rlJJfii:^ 
[0111] • ±E£Sfi&JF2fllT-f4, K«t#2 6 av</u 

s ^yzm^i-zmmzmttftmrnm 1 o 

KoT«rtKA£ft*mc0^ffl9£fr5 - fc 
10 [0112] • Sfc, M*tfK«i§a&2 0 0>»A£»O 

— /W&m&fflffll-Z*?— ^#***HUR«iiB2 Of*] 

14, _LI53£!fca. bftU«*»a<EB»^*5»*5T«ffl»« 
JE PmTAXt/^ o y h/U#ifiia«mtTA-S:±E««0»iJ» 

[0113] • ±E*Jtifi^ffi0^jBlJftL3S^*5t>Tl4. 
fartKAS**mcSr*«>5BRjc % raffi«fime*:T«IHIM 
20 K*aTniRr/Ti««BR«UEPm(OH»i:^5ll** (1 
4) ^^5«t:piafc^ HS«imcfe#^5*S 

[0114] • ±E*3IJ£«fl8-C*4, KAffiJW>«Efl« 
r^t Lt#5>*i5±ia*a (8) , (9) 

(8) . (9) S:*«Wlc«v^fc«»*(^:S^v^ 

30 [oils]- ±E#si*6?Kja-ci4, jjRwatw-aoiftB 

*Cffiw^e>*T.51BrtKAS«fl:mcS:^»J-j--<<. ±E^ 
fflfflBa t f wd Cr HI ttfl-Vltfto W: tti (ciS L fclElc^f s 

-fZZtlzX*). *^cB*J«fBUW«v *fl!>tt©«BB«»lc 

[0116] • ±E#HJ£JBl»T-f4, i3^»i 7 {^-T 

40 ^ftiii^J^O ri/2j l/n« (n :flEft) i 

[m®ofiH¥fti^0^] 

[Ull] « 1 0HJ£»ffi(CibMt5KA2g«fl:^»JSE0 



[SI 2] ^T7u — *-*<DmftW<DWrm\3o 

[13 3] JSloKjft^lBoaSBlcJiaBRAapiRSo^W 

so [[15] «HO*lk9fl|oSEBlcJ:««Wi&^W 5 
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[16] »2ojiife}Bffi«c*>jt5KA^fta-^fflygsao 

UttiMtH. 
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*n^h/U# % 2 7-±««SS^ 2 8-T««fB»; 
2 9-«M, 2 6-SMWK 3 0-ftJEttlK, 3 1- 
filgffi, 32->^Ji, 4 1-x77d^-^, 4 2 
— 9WSa-fe>1»-. 4 3-7Qy h/UtVl 4 4— * A 
A-fc^l^ 4 5-^7^«tyl 4 6-5RSU£-tr> 

iK 4 i o-«u», 4 i 2"--«Mk. 
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